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Abstract 

r 1 ■ r(ir°o ;£;;>30MeV,6»;>20°) 

We present the measurement oi the ratio R = ' ■rrn and a hrst 

measurement of the direct emission contribution for the same process. We use 328 
pb~^ of data collected at KLOE in 2001 and 2002, corresponding to about 3 million 
of -f^e3(7) events and about 9 thousand K^^^ radiative events. Our result is R= 
(924±235jai±165ysj)xlO~^ for the branching ratio and {X) = -2.3±1.3stat ±l-4sy<jt 
for the parameter describing direct emission. 

key words: direct emission 
PACS: 13.20.-v, 13.20.Eb 



1 Introduction 



The study of radiative Kl decays offers the possibility to obtain informations 
on the kaon structure and the opportunity to test theories describing hadron 
interactions and decays, Uke chiral perturbation theory (ChPT). Two differ- 
ent contribute to the photon emission, inner bremsstrahlung (IB) and direct 
emission (DE). DE is radiation from intermediate hadronic states and is thus 
sensitive to hadron structure. In /T^g DE is 1% or less of IB which diverges 
both at E*^0 and 9*^0 (photon angle w.r.t. lepton). We therefore exclude 
small angle and energy photon. To compare our result with other measure- 
ments we only retain events with E*> 30MeV and 9*> 20° [1]. We define R 
as 

„ r(/i:Js,;£);>30MeV,»;>20°) 

ms^) 
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Predictions for R ranges between 0.95x10"^ and 0.97 xlO~^ [2]. Recent mea- 
surements of R from NA48 and KTeV [3, 4] are in marginal disagreement 
between each other, so that new measurements are welcome. Following the 
authors of Ref. [2], in ChPT the structure-dependent (SD) terms are charac- 
terized by six amplitude {Vi, Ai}, which in the one- loop approximation are 
real function, almost constant over phase space. In particular, all relevant SD 
terms have a similar and simple photon energy spectrum, with a maximum 
around -E'*~ 100 MeV. This suggests to decompose the photon spectrum in 
the following manner 

'^''■-{X)f{E;) (2) 



in which the different SD contributions are summarized in the so-called dis- 
tortion function, f{E*), which represents the deviation from the pure inner 
bremmstrahlung. All the information on the structure-dependent terms is con- 
tained in the effective strength, {X), that multiplies f{E*). ChPT calculation 
at 0{p^) order from Ref [2] gives 

(X) = -1.2±0.4 (3) 



A first attemp to measure the DE contribution was performed in 2001 by 
KTeV collaboration [5], but the uncertainties due to their working hypothesis 
were too large to infer definitive conclusions on the (X) parameter. 
In our analysis we can isolate DE from IB only because we use both the energy 
spectrum and the angular distribution of the radiated photon. 



2 The KLOE detector 



The KLOE detector consists of a large cylindrical drift chamber (DC), sur- 
rounded by a fine grained lead-scintillating fiber electromagnetic calorimeter 
(EMC). A superconducting coil around the calorimeter provides an axial mag- 
netic field of about 0.52 T. 

The drift chamber [6], 4 m in diameter and 3.3 m long, is made of 58 con- 
centric rings of drift cells arranged in a stereo geometry. It is filled with a 
He, iC^Hio mixture. The spatial resolutions are axy ~0.15 mm and az —2 mm. 
The transverse momentum resolution is Cp^/px ~ 0.4%. Two-track vertices 
are reconstructed with a spatial resolution of ~ 3 mm. 

The calorimeter [7] is divided into a barrel and two endcaps. It covers about 
98% of the solid angle. It is segmented in depth in five layers, about 3Xq 
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each. The barrel is divided in 24 sectors, 5 x 12 calorimeter cells each, read 
out by photomultipliers at both ends to measure the arrival time of particles 
and to reconstruct the space cordinates. Cells close in time and space are 
grouped into calorimeter clusters. The energy and time resolutions are OeIE = 

h.l%/^E (GeV) and (Jt = 54 ps (GeV)©100 ps, respectively. The spatial 
resolution are axy ~1.3 cm and (Jz —I cm / ^ E{GeY). 

The KLOE trigger [8] uses calorimeter and chamber information. For this 
analysis, only the calorimeter signals are used. Two energy deposits above 
threshold {E > 50 MeV for the barrel and E > 150 MeV for the endcaps) are 
required. Recognition and rejection of cosmic-ray events is also performed at 
the trigger level. Events with two energy deposits above a 30 MeV threshold 
in the outermost calorimeter plane are rejected. The 328 pb~^ of data used in 
this analysis, taken in 2001 and 2002, are divided in 14 periods of about 25 
pb~^/period. For each data period we have a corresponding period simulated 
with Monte Carlo (MC) with about the same statistic. 



3 Monte Carlo Generator 



In the KLOE MC only radiation from inner bremsstrahlung is described, so we 
need a Monte Carlo to describe the photon spectrum from direct emission. For 
this purpose we use a Monte Carlo generator {Kubis generator) based on the 
code provided by one of the autors of Ref. [2], which implements their 0{p^) 
calculation. We generate DE events folded with the KLOE reconstruction MC. 
The accuracy of KLOE Monte Carlo in describing the photon spectrum from 
IB is at level of ~l-2% (an appropriate level for many KLOE measurements). 
In particular, the KLOE MC generator avoids the problem of the infinite value 
for the total decay width for a single photon emission by re-summing, in the 
limit of soft photon energy, the probabilities for multiple photon emission to 
all order in a [9] . Unfortunately, this accuracy level is of the same order of DE 
contribution, which is ~1% of IB one. From a point of view of the measurement 
of R, this could introduce only a ~1% error. On the other side, a fit counting 
procedure based on a 1%-biased IB distribution could introduce up to ~100% 
error counting in K^^^ from DE events. For this reason in this analysis we use 
the Kubis Monte Carlo generator also to describe the photon spectrum from 
IB. 
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4 Analysis 



The criteria used to select -f^'e3(^) events, briefly summarized below, are the 
same described in Ref. [10]. Candidate Kl events are tagged by the pres- 
ence of a Ks — > Ti'^Tx'^ decay Fig. 1 shows that the tagging efliciency, about 
66%, is almost independent of the photon energy. The Kl is searched along 
the direction of its momentum {tagging line), reconstructed from the decay 
Ks — > TT+TT". -f^e3(^) evcuts are then selected using appropriate kinematical 
variables of the decay and electron identification by time of flight (TOF). 
We have about ~3 million of K2nu.\ events with a contamination of 7x10^'^, 




20 40 60 80 100 120 140 160 180 

E • (MeV) 



Fig. 1. Tagging efficiency from Monte Carlo simulation as a function of plioton 
energy. Full statistic 

mainly due to -f^°3(^) events. 

To select K^^„^ events we search for a photon cluster, i.e., a cluster in the 
calorimeter not assigned to any track. The arrival time of each photon gives 
an independent determination of the position of the K^ vertex, X^r, the so- 
called neutral vertex (NV). The method is fully described in Refs. [7] and [11]. 
The position of the K^ vertex is assumed to be along the K^ line of flight. 
We require that the distance, d^c^ between the position of the neutral 
vertex and the position Xc of the Ki vertex determined with tracks, to be 
within 8cr. In case of more than one photon candidate, we choose the closest 
to the Kl charged vertex. To evaluate the photon energy we use the charged 
track momenta and the photon cluster position, X^u. By solving equation 4 
below in the hypothesis of neutrino zero mass, we extract the photon energy 
with a resolution of ~ 1 MeV. This resolution is a factor ~ 10 better than that 
obtained using the energy deposit information of the calorimeter. 



Pu=PK-p-K-Pe- Pj 



N 



IX, 



clu 



X 



(4) 



N\ 
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In this equation p^, px, Pn, Pe and p^ are the neutrino, kaon, pion, electron 
and photon momentum, respectively. 

Fig. 2(a) gives the selection efficiency for the signal as taken from Monte Carlo 
simulation. The main background contribution comes from K^^ events with an 
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(a) (b) 

Fig. 2. Monte Carlo simulation: (a) K^^^ signal efficiency (one run period); (b) 
applied cut to remove accidentals. 

undetected soft-photon to which a cluster from accidentals has been wrongly 
associated. This background is strongly reduced by requiring >25 MeV 
and ^ciu -^^"^ < ^ciu - 15 MeV where E^^ is the energy of the photon cluster 
and E}^^ is the reconstructed photon energy in the laboratory system. This 
cut is shown in Fig. 2(b). We obtain a factor 10 in background reduction with 
~6% loss in signal efficiency. 

Background from Kl tt^tt^tt^ and Kl — ^ vr/iz/ events after signal selection 
is at level of ~ 4.2% and ~ 2.5%, respectively. As shown in Fig. 4, this back- 
ground overlaps with the signal of interest DE, so we need to remove it as much 
as possible. To remove both tt'^tt'tt'^ and Kl — > tt/xz/ we use a neural 

network (NN). To remove Kl tt^tt^tt^ events we use a NN based on the 
photon energy and angle {w.r.t. the lepton), the track momenta, the missing 
momentum and M^^, the invariant mass of photon- neutrino pair. To remove 
Kl 7r//z/ events we use a NN based on the track momenta, the calorimet- 
ric energy deposit and the cluster centroid position. Appropriate cuts on the 
NN output give a background reduction from 4.2% — > 0.4% and from 2.5% 
1.4%, respectively for Kl 7r+7r^7r° and Kl ~^ irfiu, with a signal loss of 
10%. To check the data-MC agreement, to calibrate the MC position and 
correct the photon selection efficiency we use Kl — > vr+7r~7r° decay events as 
a control sample. These events are selected using a tight kinematical cut in 
the variable -Emiss ~ Pmiss ~ ''^l-o in the hypothesis of two pion tracks. Further, 
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we require the presence of a cluster (E>60 MeV) not associated to any track, 
corresponding to one of the two photons from 7r° decay. This high energy pho- 
ton is used to tag the presence of the second photon. We select about 350,000 
Kl — > 7r"^7r~7r'^ events with a purity of 99.8%. 

As a first check of Data-Monte Carlo agreement we compare the energy res- 
olution of the photons. This can be done because in this control sample we 
estimate the energy of the second photon (the tagged photon) exactly in the 
same way (Eq. 5) than in our K^^^^ signal selection: there we do not detect the 
neutrino, here we ignore the hard photon (the tagging photon). Squaring the 
equation below the (second) photon energy is extracted. 

P^-hard = PK -Ptt -Ptt - P-y ] P7 = -^7 ^ (5) 

l-^clu — Xn\ 




Fig. 3. From Kl vr+vr vr*^ control sample: (a) photon energy resolution; (b) d]\fc 
distribution in the central part of the drift chamber after correction. 



The photon energy resolution is evaluated with respect to a value computed 
in a more accurate method using the complete hard photon informations and 
closing in this way the kinematic. In Fig. 3(b) the residual is shown for data 
and Monte Carlo and a good agreement results. 

Further, we use the control sample from Kl tx^ti'ti^ to evaluate the dis- 
tance (Inc and its resolution 0"^^^, in order to correct (Inc and cxd^^ in MC 
simulation. Because of the use of -Eciu to remove accidentals, we also use this 
control sample to check the calorimeter energy response. The MC energy re- 
sponse is about ~2 MeV lower than data. To a good approximation, this bias 
is independent of the energy. 

Finally, we evaluate the photon selection efficiency from data and MC in this 
control sample and use their ratio to correct photon selection efficiency in MC 
simulation. The correction is of the order of a few percent. 
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5 Fit 



In order to count K^^^ signal events we fit the Monte Carlo spectra fi{E*, 9*) 
to tlie data {i = 1,2,3,4 respectively for IB signal, DE signal, K^^^ont- 
of-acceptance {E* < 30MeV or 9* < 20°) and physical background from 
Kl — >■ 7r"'"7r~-7r° and Kl — > irfiu events). The four distributions used as inputs 
in the fit are shown in Fig. 4. Actually it is possible to measure R by using the 
energy spectrum of the photon only (we do it and the result matches), but in 
this case there is no sensitivity to the presence of a DE term. Only a simultane- 
ous use of the energy E* and the angle 9* can disentangle the small DE signal. 
To check the fit stability as a function of run period we fix the background 




Fig. 4. Monte Carlo distribution, 9* (deg) vs E* (MeV). From left-top: IB sig- 
nal, DE signal, -R'g3^-out-of-acceptance and background from Kl — tt^tt~tt^ and 
Kl — > TTfii'. The vertical scale is in arbitrary units. 

component, otherwise the fit could not converge (too low background statis- 
tic) and we do not use the DE shape (no sensitivity in a single run period). 
The stability is good (x^/dof = 9/13). Then we fit all data simultaneously: 
free parameters of the fit are the population for IB signal, DE signal, K^^^- 
out-of-acceptance. We fix the background contribution of Kl iT~^7r~7T^ and 
Kl nfiu from MC. As a check we also perform the fit letting all parameters 
free. The result matches well but there is a loss in statistical accuracy. For this 
reason the background fraction is not a free parameter of the fit. The result of 
the fit and the residual are shown in Fig. 5. The two-dimensional MC input 
shapes are arranged as 8 ^ — slices one dimensional hystograms. Each slice 
covers 20 degrees, from 20° to 180°. The counting result and the correlations 
of parameters are given, respectively, in Table 1 and Table 2. The is good. 
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(a) (b) 

Fig. 5. (a) Fit: dots are data, dark grey is the signal (IB+DE), light grey is the 
Kgg^-out-of-acceptance; (b) Fit residual. 

giving x^/dof=60/69). We get the -f^e3(7) events by counting events after TOF 





Counting result 




IB signal 


i^'g3^-out-acc 


DE signal 


counting 


9083 


6726 


-102 


error 


213 


194 


59 



Table 1 

Conting from the fit with the statistical error 





Correlation coefficients 


Par 


1 


2 


3 


1 


1.000 


-0.586 


-0.254 


2 




1.000 


-0.022 


3 






1.000 



Table 2 

Correlation coefficients: 1: IB signal, 2: K^^^ -out-of- acceptance, 3: DE signal 

selection, bakground subtraction and efficiency correction. After getting K^^^ 
signal from ffi counting, taking into account the efficieny and all the efficiency 
corrections we measure: 



R 



BR{K%^- E; > 30 MeV, 9; > 20°) 



BR(i^o 



(924 ± 23) X 10" 



(6) 



e3(7)^ 



Although the for the ffi is not bad if we use the KLOE MC to describe the 
IB spectrum (x^ probability ~55%) the Kubis Monte Carlo is slightly better, 
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giving a 77% probability. 



6 Systematic uncertainties 

We estimate all the systematics uncertainties by varying the cuts. Tracking, 
clustering, track-to-cluster association, NV acceptance and analysis cuts all 
depends on some parameters which define our signal. Any variation on these 
parameters produces a variation on the result. In the following we list the 
absolute variation of 10^ x R. 

Tagging. We tag the Kl requiring that Ks alone satisfies the calorimeter trig- 
ger with the presence of two clusters from Ks —>■ 7T~^7r~ associated with fired 
trigger sector (autotrigger). We observe a change of 4. 

Tracking. The most effective variable in the tracking candidates definition 
is dc, the distance of closest approach of the track to the tagging line. We 
vary dc by a factor of two. We re-evaluate for each different configuration the 
tracking-efficiency correction, which is run-period dependent. The uncertainty 
on the tracking efficiency correction is dominated by sample statistics. We 
observe a change of 1.5. 

Clustering. The most effective variable in the definition of TCA association is 
the transverse distance, dtc±. We vary the cut on dtc± from 15 cm to 30 cm, 
corresponding to a change in efficiency of about 17%. We re-evaluate for each 
different configuration the clustering efficiency correction, which is run-period 
dependent. Also in this case the uncertainty on the clustering efficiency cor- 
rections is dominated by sample statistics. We observe a change of 5.5 in the 
result . 

Kinematic cuts. We apply loose kinematic cuts. When varying this cut negli- 
gible variation for the results are found. 

To f cut. Inclusive -f^e3(7) sample is identified also by using time of fiight (TOF). 
We use a 2-a cut. After varying this cut by 30% we observe a change of 1.3. 
Momentum mis- calibration and resolution. The effect of the momentum scale 
and the momentum resolution have also been considered. We conservatively 
assume a momentum scale uncertainty of 0.1% We observe a change of 3 for 
R. 

We also investigate the effect of momentum resolution by changing its value 
of ± 3%, corresponding to a worst x^- The variation on the result is 7.2. 
Fiducial volume. We reduce the fiducial volume by a 20%. This produces a 
variation of 3. 

Rejection of accidentals. We vary the sliding cut used to remove wrong associ- 
ations of accidental cluster (see Fig.2). Varying this cut we have a background 
variation of a factor of two. The change in the result is 5.2. 
NV acceptance. We search neutral vertex within a well defined sphere centered 
around Xq- We vary the dimension of its radius by a factor of two. We observe 
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a change of 2.9 in the result. 

Background. In a very conservatively way we remove the cut on NN output: in 
this way the background level increses almost of a factor of four. We observe 
a change of 9 for the result. 
All systematic errors are summarized in Table 3. 



Source 


10^ X AR 


A{X) 


Tagffine: 


4.0 


0.7 


Tracking 


1.5 


0.8 


Clustering 


5.5 


0.1 


Kinematic cuts 


~0 


~0 


TOF-cut 


1.3 


0.5 


p-miscalibration 


3.5 


0.2 


p-resolution 


7.2 


0.4 


Fiducial volume 


3.0 


0.5 


Rejection acc. 


5.2 


0.4 


NV acceptance 


2.9 


0.3 


BKG Reject .by NN 


9.0 


0.1 


Total 


15.5 


1.4 



Table 3 

Summary of the absolute systematic uncertainties on R and {X) 



7 Results 

Our final result for R is 

i?= (924±23stat±16syst) X 10-^ (7) 

We also estimate the (X) parameter, defined in the equation 2. Starting from 
counting result for IB and DE, taking into account the difference in the effi- 
ciency (IB efficiency is about ~20% higher than DE efficiency), including all 
systematics we measure 

{^)m.eas = —2.3 ± 1.3stat ± l-4syst (8) 



in agreement with evaluation. The systematics on (X) are evaluated as 

for R. The different contributions are listed in Table 3. The presence of DE 
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contribution reduces the value of R of about 1%. The correlation between R 
and (X), including also systematics, is 3.9%. 



2 




-2 



-6 
-8 

900 950 

Fig. 6. 1-a confidence level for R and (X) measurement. 



8 Conclusion 

Two different components contribute to the photon emission in -^^£37 events, 
the inner bremsstrahlung and the direct emission. The latter describes photon 
radiation from intermediate hadronic states, giving in this way new informa- 
tion on the hadronic structure of kaons. Using E* and 6* variables KLOE 
measured the width for Kl 7r^e'^u{'y) for E* > 30MeV and 9* > 20° to 
the width for inclusive K'^^^y The DE emission contribution originates from 
the interference with IB, resulting in a negative effective strength, {X). KLOE 
measurement of {X) is the the first attempt to measure direct emission con- 
tribution in K^2^ process. At this stage, the KLOE measurement of R (3% 
accuracy) is not sufficient to solve the experimental disagreement between 
NA48 and KTeV measurement. 
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